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Mathieu ordinary differential equation is of Fuchsian types with the two regular and 
one irregular singularities. In contrast, Heun equation of Fuchsian types has the three 
regular and one irregular singularities. Heun equation has the four kind of confluent 
forms: (1) Confluent Heun (two regular and one irregular singularities), (2) Doubly 
confluent Heun (two irregular singularities), (3) Biconfluent Heun (one regular and 
one irregular singularities), (4) Triconfluent Heun equations (one regular and one 
irregular singularities). Mathieu equation in algebraic forms is also derived from the 
Confluent Heun equation by changing all coefficients 6 = 'j = ^, e = 0, a = q and 
q = [see (67) in |arXiv:1303"M76] 

In this paper I will apply three term recurrence formula by Choun [a rXiv: 1303.0806] 
to the power series expansion in closed forms of Mathieu function (infinite series) 
and its integral forms including all higher terms of A^s. One interesting observation 
resulting from the calculations is the fact that a Modified Bessel function recurs in 
each of sub-integral forms: the first sub-integral form contains zero term of A'^s, the 
second one contains one term of A^s, the third one contains two terms of AnS, etc. 
Section VI contains two additional examples of Mathieu function. 

This paper is 5th out of 10 in series "Special functions and three term recurrence 
formula (3TRF)". See section VI for all the papers in the series. Previous paper in 
series deals with asymptotic behavior of Heun function and its integral formalism 



by Choun (arXiv: 1303.0 876J . The next paper in the series describes the power series 
expansion in closed forms of Lame function and its integral forms in the algebraic 



form by Choun [arXiv: 1303.0873] . 
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I. INTRODUCTION 

One example of three therm recursion relations is the Mathieu equation, introduced by 
Mathieu (1868)^, while investigating the vibrating elliptical drumhead. Mathieu equation, 
known for the elliptic cylinder equation appears in diverse areas such as astronomy and 
physical problems involving Schrodinger equation for a periodic potentials^, the parametric 
Resonance in the reheating process of universe^^, and wave equations in general relativity^^, 
etc. Mathieu function has been used in various areas in modern physics and mathematics.-"— 

Unfortunately, even though Mathieu equation has been observed in various areas men- 
tioned above, there are no power series expansion in closed forms and its integral formalism, 
analytically. Mathieu function has only been described in numerical approximations (Whit- 
taker 1914i^, Frenkel and Portugal 2001^^). Sips 1949^^, Frenkel and Portugal 2001-^ argued 
that it is not possible to represent analytically the Mathieu function in a simple and handy 
way. 

I will construct the power series expansion of Mathieu function in closed forms and its 
integral formalism analytically using three-term recurrence formula^ and the same method 
I used in analyzing Heun function.— 

Mathieu equation is 



where A and q are parameters. This is an equation with periodic- function coefficient. Math- 
ieu equation also can be described in algebraic forms putting x = cos^z: 



This equation has two regular singularities: x=0 and x=l; the other singularity x=oo is 
irregular. Assume that its solution is 




(1) 





oo 




(3) 



n=0 



Plug dS]) into 02]). 




■n> 1 



(4) 



where. 



A. 



4{n + vf - (A + 2g) 



(5a) 



■n 



2(n+l + z/)(2(n + i^) + l) 
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4g 

^" " 2(n + l + z/)(2(n + z/) + l) ^^^^ 



l-(A + 2g) 

ci = Co = ^ Co (5c) 



We have two indicial roots which are z/i = and 1^2 = \- As we see fl5bp . there is no way to 
make Bn term terminated at certain value of n. Because the numerator of (15bp is just consist 
of constant q parameter. So there are only two kind of power series expansions: infinite series 
and the polynomial which makes term terminated. In this paper I construct an analytic 
solution of Mathieu equation for the infinite series. And in the future I will construct its 
analytic solution for the polynomial which makes An term terminated. 



II. POWER SERIES EXPANSION FOR INFINITE SERIES 

In Ref. Q, the general expression of power series of y{x) for infinite series is 

yi.x) = ^yn{x) = yo{x) + yi{x) + y2{x) + y^ii^x) ^ 

{00 Ao— 1 \ 00 ^' io— 1 00 / 42 — 1 \ ^ 

E n ^-^+1 + E n B2n^^ E n ^-3+2 

«0=0 \ii=0 / 40=0 L n=0 12=10 \i3=io / ) 

00 f 00 iQ — l TV— 1 / 00 *2fc~l 

+ E ] E 1 "^2*0 Yl B2h + l W ( E "^2i2fc+fc W 52i2,+,+(fc+l) 
JV=2 I io=0 L ii=0 k=l \«2fc=«2(fc-i) «2fc+l=«2(fc-i) 

X E I n 5-..,,+(TV+i)) (6) 



«2JV— «2(]V-1) \ «2JV+1— «2(JV-1) 



Substitute (I5aj) - fl5c|) into ([6]). In this article Pochhammer symbol (x)„ is used to represent 
the rising factorial: (x)„ = ^^f^^- And the general expression of power series of y{x) for 
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infinite series is 



oo 

y(^) = ^Vni^) = yo{x) +1/1 (x) +y2ix) +?/3(a;) H 

n=0 




2 

>p f >p f (^o + |)^-^(A + 2g) 1 
^ Ifc^o I (^0 + i + |)(^o + i + f ) (1 + f ).o(! + |).o 



X 



I I , X- fa + I + 1)^ - ^(A + 2g) (1 + I + |).,.,(| + I + 

tiU.i:.(^^ + t + i + t)^^^ + t + i + f) (i + l + fk(! + l + fk 



^ri — 1n — l 



(7) 



As z/=0 and co= ^pAr- in ([7]), 



= I X = cos'^z, 7] = -qx^ 




fa + |)fa + i) 





n=2 L jo=0 
n— 1 f oo 

n E 



^(A + 2g) 1 

fa + |)fa+i) (l)io(!).o 

fa + |)^-^(A + 2g) + + 
fa + | + i)fa + | + i) (l + |k(| + i.. 



(1 + 2 (4 2 




^^"-^^l-^if^^^-^K^ (8) 



([8]) is the 1^^ kind of independent solution of Matliieu function for infinite series. 



5 



Asi/ = iand co= (ig)Vr(|) in©, 
y{x) = Ms (x = cos^z,r] = -^QX^ 



'7\ /3^ 



U^o(t)-W- t^ol (^o + !)(^o + i) (|k(ik.^a(|k(|k 

, ff^ [ (^o + i)^-^(A + 2g) 1 

"tt7 fa + i + r-^(A + 2g) (! + |k-.(l + |k-. 

^ (^. + | + |)(^. + | + f) (! + fk(i + fk 

k-.(l + t).. 



X 



X. (4 + 2)^n-l(l + 1 I 

^ji — ^n — l 



is the 2'^'-^ kind of independent solution of Mathieu function for infinite series. 



III. INTEGRAL FORMALISM FOR INFINITE SERIES 



There is a generahzed hypergeometric function which is 



2 ' 2/b-iV4 ' 2 ' 2/«j-i 



(i + l + lk(! + l + ^k 



ij=ij^i ' 2 ' 2^*jV4 ' 2 ' 2^*J 

00 ^ 



By using integral form of beta function, 



+ | + + = l\,1^-^-^^^^^{l-tjy (lla) 



B{^,-l + '- + -,l + l] = I du.ul' '-^-^{l-u,r (lib) 







Substitute (fTTal) and flllbl) into ([TO]). And divide (ij-i - ^ + f + f)(^i-i + i + f ) into Gj. 

J- \ ^ '^-^ ^ 2 ^ 2^j-lU ^ 2 ^ 2Aj-i i. 



(^.-i-i + i + f)(vi + i + f).,ft, (i + i + fk(! + l + 



6 



The modified Bessel function is 



~ 1 /x\2;+" (^) 



replaced a and x by and 2^Jr]{l — tj){l — Uj) in ([T3 



(14) 

Substitute (IHD into flT2D. 



(-*- + 2 + 2)^-1(4 + 2 + 2)^-1 i,- 

(^.-i-i + i + f)(^.-i + i + f).X. (1 + 1 + 1)^1 + l + fk 



^0 



Defined by 

b 

Wa,b = { l=a (16) 

T] only if a > 6 

Substitute f[T^ into ([7]) by applying f lTB]) . The integral formalism of the function ?/(x) for 
infinite series is 



00 ^ n— 1 I' „x 

+E Hi / ^^-'^^n 

n=l I fc=0 I 



\ + \(n-k+v) 1^ J _-l+^{n-k+u) 



-i(n-l-fc+i/) / x2 i{n-l-fc+!.) 1 , . 

^n4n [Wn-k,ndy,„.,J <_fc,„ -— (A + 2g) 
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Let z/ = and Cq = ^prrr in ( !T7|) . 
y{x) = Mf (x = cos^z, 7] = -^qx^ 



_n l-'-JioU^'w) «-i i,-n >^0 ^0 



Iq {^\J Wn+i^k,niX -tn-fc)(l -Mn-fe)^ (18) 

1 ]^ 1^1 



^n—k,n 



replaced a and x by — ^ and in (IT^ . 

4 °^ 1 „-i ^1 



— - ^ 1 — — 

- E 171);"' - fi|)m) - '''^ 



Similarly, 



1 1 

"8 00 ^ ,,,8 /•! 



8 1 8 /* 1 



Substitute (1191) and (EOl) into fITSll. 



= I X = cos'^z, 7] = -qx" 



oo f n— 1 f „i „x 

x/o |^2y^u;„+i_fe,„(l -t„_fe)(l - Un-k)^ (21) 



X <j ^„-M ' '^^ {Wn-k,nd^,,,_,^S <-k,l ''^ - ^(^ + 2g) 1^ (2^^) 



( 12T|) is the integral formalism of the first kind of independent solution of Mathieu function 
for infinite series. 

Let z/ = i and co = ^ (^g)^ in (HZ])- 

y{x) = Ms = cos^z, r] = -qx'^ 

1 ^ oo -I oo ( n— 1 f pi nl 



,„ _n ll^ol-'-Jjo „_i i,_n -^0 Jo 



x/o ^2y^w„+i_fc^„(l - t„_fc)(l - Un-k)^ (22) 

X l^n-tn '' (^n-fc.n-^^n-fc.nj^^-Tn''^ " ^(^ + 2g)| ^^-(Jj-JT^'^l,n ^ 



*0 11^" 

20 
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replaced a and x by | and 2^7/ in (fT3!) . 



r(!)t^^!(!)/ r(i)r(|) 

Similarly, 



1 1 

oo 



8 1 8 /• 1 



Substitute (El and (El into (l22D. 



y{x) = Ms ( X = cos^z, ?7 = -ga;' 



{oo f n—1 f pi pi 
n=l I fe=0 I ^0 

x/o t(;„+i_fc,„(l - t„-_fe)(l - Un-k)^ (25) 



( 125|) is the integral formalism of the second kind of independent solution of Mathieu function 
for infinite series. As we see ( 12T1) and ( l25ll . Integral forms of Mathieu function are consist 
of modified Bessel functions surprisingly. We can transform Mathieu function from these 
integral forms to other well-known special functions: Kummer function, Legendre function, 
Hypergeometric function, Laguerre function etc. 

For the special case, under the condition A ^ g > 0, the approximate solutions of the 
Mathieu equation in (jH]) and (E]) are 

lim Mf \ X = cos^2, T] = -qx' 
A»g>0 M ' ' 4^ 



r(!) 1,^0 1(^0 + i)(^o + i)(ik(!k .l^„(|k(!k' IV 4^" 

1 



oo f oo 



(io + 5)(«o+i)(lk(|)i 



X 



X y (26a) 
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also, 



lim M, ( X = cos^z, 7? = -qx'^ 

X»q>0 M '4^ 




°° fl). (3.). 




1 1 



X 



(^o + !)(^o + i)(|k(lk 

TT j 1 (f + |)»fc-i(l + |;tfc- i 

oo 5 I n\ f 1 _|_ 

4 2 V ^ 2 



— ^n— 1 




And integral forms of approximated Mathieu function of f l26al) and f l26bl) are 



lim Mf X = cos z,n = -qx^ 
A»g>o ■'V 4 



, -l+|(n-fc) 



^ oo f n-i f „i 

"f(i)f(i)g|n|x 

X / dv,{i-viru-^^^-\\[-^^x 



(27a) 



and. 



lim Mo X = cos z,ii = -qx' 
A»g>0 'V '4^ 



1 



r(i)r(^ 



oo ( n — 1 

En 



^ ,-l+|(n-fe) /"^ _3 + i(„_fc) 



^4/ n=l I, fc=0 



-1 



A 



(27b) 
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IV. ASYMPTOTIC BEHAVIOR OF THE FUNCTION y{x) AND THE 
BOUNDARY CONDITION FOR x = cos^ z FOR THE CASE OF INFINITE 
SERIES 



As n 1, ( I5all and ( l5b1l are 



\imAn = A = l (28a) 

n>l 

And, 

lim 5„ = 5 = ^7 (28b) 

n>l n"' 

As n ^ 1, (]28bp is extremely smaller than fl28ap as it's compared to each other. Then, put 
fl28al) with 3^ = into i^. 

Cn+l ^ Cn (29) 
oo 

Plug (!29l) into the power series expansion where ^^c^x", putting Cq = 1 for simplicity. 

n=0 

oo ^ 

lim ?/(x) = = where < a; < 1 (30) 

"'^ n=0 

As we see ( 130|) . x = cos^ 2; can not be 1 for the case of infinite series. If it does't, the function 
y{x) will blows up as n ^ 1 and x = 1. 

I show the power series expansion in closed forms of Mathieu function for the infinite 
series in this paper analytically. Also, I derive integral forms of Mathieu function from 
power series expansion. It is quiet important that integral forms of this function, because 
we can investigate how this function is associated with other well known special functions; 
Bessel function, Kummer function, hypergeometric function, Laguerre function and etc. And 
we can obtain orthogonal relations, normalized physical factors and expectation values of 
any physical quantities from the integral forms of its function analytically. 



V. APPLICATION 

1. By using my methods on the above, we can apply the power series expansion of the 
Mathieu function and its integral forms into various modern physics areas.-""— For example, 
in general Mathieu equation arises from two-dimensional vibrational problems in elliptical 
coordinates with physical points of a view-^. Its equation is derived from the Helmholtz 
equation in elliptic cylinder coordinates by using the method of separation of variables (see 
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(5)-(7) in p. 610 in Ref. llTl ). Using the power series expansion of Mathieu equation, it might 
be possible to obtain specific eigenvalues for the wave equation in vibrational systems. Using 
the integral forms of Mathieu equation, it might be also possible to construct the normalized 
wave functions for the entire region. 

2. In "Examples of Heun and Mathieu functions as solutions of wave equations in curved 
spaces"—, "Dirac equation in the background of the Nutku helicoid metric"—, two authors 
consider the Dirac equation in the background of the Nutku helicoid metric in five dimen- 
sions. They obtain the solution for the four different components by using the Newman- 
Penrose formalisn>22(see (8a)-(8d) Ref. [lol). And they separate uncoupled equations for the 
lower components into two ordinary differential equations: angular equation is of the Math- 
ieu type and radial equation is of the double confluent form which can be reduced to the 



Mathieu equation (see (11), (13) Ref. llQl ). Using power series expansions and integral forms 



of Mathieu equation, it might be possible to obtain eigenvalues and normalized wave func- 



tions at various regions. In Ref. |28|, I will investigating power series expansions in closed 
forms and integral forms of Lame function in the algebraic form and its asymptotic behaviors 
analytically. 



VI. SERIES "SPECIAL FUNCTIONS AND THREE TERM 
RECURRENCE FORMULA (3TRF)" 

This paper is 5th out of 10. 

1. "Approximative solution of the spin free Hamiltonian involving only scalar potential 
for the q — q system" - In order to solve the spin- free Hamiltonian with light quark masses 
we are led to develop a totally new kind of special function theory in mathematics that 
generalize all existing theories of confluent hypergeometric types. We call it the Grand 
Confluent Hypergeometric Function. Our new solution produces previously unknown extra 
hidden quantum numbers relevant for description of supersymmetry and for generating new 
mass formulas. 

2. "Generalization of the three-term recurrence formula and its applications"^^ - Gener- 
alize three term recurrence formula in linear differential equation. Obtain the exact solution 
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of the three term recurrence for polynomials and infinite series. 

3. "The analytic solution for the power series expansion of Heun function" - Apply 
three term recurrence formula to the power series expansion in closed forms of Heun function 
(infinite series and polynomials) including all higher terms of AnS. 

4. "Asymptotic behavior of Heun function and its integral formalism" ,— - Apply three 
term recurrence formula, derive the integral formalism, and analyze the asymptotic behavior 
of Heun function (including all higher terms of y4„s). 

5. "The power series expansion of Mathieu function and its integral formalism" ,— - Apply 
three term recurrence formula, analyze the power series expansion of Mathieu function and 
its integral forms. 

6. "Lame equation in the algebraic form"— - Applying three term recurrence formula, 
analyze the power series expansion of Lame function in the algebraic form and its integral 
forms. 

7. "Power series and integral forms of Lame equation in the Weierstrass's form and its 
asymptotic behaviors"— - Applying three term recurrence formula, derive the power series 
expansion of Lame function in the Weierstrass's form and its integral forms. 

8. "The generating functions of Lame equation in the Weierstrass's form"=2fi - Derive the 
generating functions of Lame function in the Weierstrass's form (including all higher terms 
of A^s). Apply integral forms of Lame functions in the Weierstrass's form. 

9. "Analytic solution for grand confluent hypergeometric function"— - Apply three term 
recurrence formula, and formulate the exact analytic solution of grand confluent hyperge- 
ometric function (including all higher terms of A^s). Replacing /i and eu by 1 and —q, 
transforms the grand confluent hypergeometric function into Biconfluent Heun function. 

10. "The integral formalism and the generating function of grand confluent hypergeomet- 
ric function" - Apply three term recurrence formula, and construct an integral formalism 
and a generating function of grand confluent hypergeometric function (including all higher 
terms of Ans). 
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